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ABSTRACT 



The low energy (< 10 keV) X-ray emission of the Soft Gamma-ray Repeater 
SGR 1806-20 has been studied by means of four XMM-Newton observations car- 
ried out in the last two years, the latest performed in response to a strong se- 
quence of hard X-ray bursts observed on 2004 October 5. The source was caught 
in different states of activity: over the 2003-2004 period the 2-10 keV flux dou- 
bled with respect to the historical level observed previously. The long term raise 
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in luminosity was accompanied by a gradual hardening of the spectrum, with the 
power law photon index decreasing from 2.2 to 1.5, and by a growth of the burst- 
ing activity. The pulse period measurements obtained in the four observations 
are consistent with an average spin-down rate of 5.5 xlO -10 s s _1 , higher than the 
values observed in the previous years. The long-term behavior of SGR 1806-20 
exhibits the correlation between spectral hardness and spin-down rate previously 
found only by comparing the properties of different sources (both SGRs and 
Anomalous X-ray Pulsars). The best quality spectrum (obtained on 6 Septem- 
ber 2004) cannot be fitted by a single power law, but it requires an additional 
blackbody component (kT BB =0.79 keV, R B b = 1-9 (d/15 kpc) 2 km), similar to 
the spectra observed in other SGRs and in Anomalous X-ray Pulsars. No spectral 
lines were found in the persistent emission, with equivalent width upper limits in 
the range 30-110 eV. Marginal evidence for an absorption feature at 4.2 keV is 
present in the cumulative spectrum of 69 bursts detected in September-October 
2004. 

Subject headings: stars: individual (SGR 1806-20) — stars: neutron — X-rays: 
bursts 



1. Introduction 

The high-energy sources known as Soft Gamma-ray Repeaters (SGRs) are probably one 
of the most intriguing manifestations of young neutron stars. They were first discovered 
as transient phenomena through the observation of short (< 1 s) gamma-ray bursts. The 
detection of several bursts coming from the same sky directions, coupled with their softer 
spectra, clearly set them apart from the standard gamma-ray bursts and led to the definition 
of this small class of sources. Only three confirmed SGRs are known in our Galaxy and one 
in the Large Magellanic Cloud (see Hurley 2000 for a review). 

The nature of SGRs bursts remained a mystery for many years. The 8 s periodicity 
seen during an exceptionally bright flare from the LMC SGR 0525-66 (Mazets et al. 1979), 
as well as the spatial coincidence of this SGR with the supernova remnant N49 (Cline et 
al. 1982), suggested that neutron stars could be involved in SGRs, but a confirmation had 
to wait for the discovery of their persistent counterparts in the classical X-ray range (<10 



1 Bascd on observations obtained with XMM-Newton, an ESA science mission with instruments and 
contributions directly funded by ESA Member States and NASA 
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keV). At these energies SGRs are generally observed as pulsating sources with periods of 
several seconds, secular spin-down in the range ~10~ n -10 -10 s s _1 , and X-ray luminosity 
orders of magnitude larger than their rotational energy loss (see, e.g., Woods & Thompson 
2004 for a recent review). 

The SGRs properties are successfully explained by the "magnetar" model (Duncan & 
Thompson 1992, Thompson & Duncan 1995). Magnetars are neutron stars in which the 
dominant source of free energy is a very intense magnetic field (B~10 14 -10 15 G), rather 
than rotation as in ordinary radio pulsars. Such a high magnetic field can be produced by 
an efficient dynamo mechanism that operates if the neutron star is born with a very short 
(~1 ms) period (Thompson & Duncan 1993). In the magnetar model the short bursts are 
produced by small cracks in the neutron star crust, driven by the magnetic field diffusion 
(Thompson & Duncan 1995), or, alternatively, by the sudden loss of magnetic equilibrium, 
through the development of a tearing instability (Lyutikov 2002, 2003). The much more 
energetic flares sporadically seen in SGRs are the result of global reconfigurations of the 
neutron star magnetosphere. 

In this work we concentrate on SGR 1806-20, the most prolific SGR, which showed 
several periods of bursting activity since the time of its discovery in 1979 (Laros et al. 1986). 
Its persistent X-ray counterpart was discovered with the ASCA satellite (Murakami et al. 
1994) and studied in more detail with RXTE and BeppoSAX. The RXTE observations 
led to the discovery of pulsations (period P=7.47 s and period derivative P=8xl0~ n s 
s -1 ; Kouveliotou et al. 1998) and were subsequently used to study in detail the timing 
properties of the source, such as the long term P variations (Woods et al. 2000, 2002) and 
the evolution of the pulse profile (G6gu§ et al. 2002). The best X-ray spectra reported 
to date were obtained by BeppoSAX in October 1998 and March 1999 (Mereghetti et al. 
2000). These showed a spectrum equally well described in the 2-10 keV range by a power 
law with photon index T=1.95 or by a thermal bremsstrahlung with temperature kT t6 =ll 
keV. Similar flux values were measured by all the above satellites, indicating a fairly stable 
luminosity of ~3xl0 35 (d/15 kpc) 2 erg s^ 1 (2-10 keV) in the period 1993-2001. Thanks to 
the unprecedented imaging capabilities of the INTEGRAL satellite in the hard X-ray/soft 
7-ray range, the study of the persistent emission from SGRs has been recently extended to 
higher energies: observations of SGR 1806-20 carried out in 2003-2004 showed a power law 
spectrum with photon index T ~1.5 1.9 extending up to 150 keV (Mereghetti et al. 2005a; 
Molkov et al. 2005). 

During the last two years SGR 1806-20 displayed a gradual increase in the level of 
activity, as testified by the rate at which bursts were emitted and by an increase of the soft 
and hard X-ray luminosity (Woods et al. 2004, Mereghetti et al. 2005a), which culminated on 
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2004 December 27 with the emission of the first giant flare seen from this source (Borkowski 
et al. 2004, Mazets et al. 2004). 

Here we report on a series of XMM-Newton observations of SGR 1806-20, obtained 
from April 2003 to October 2004, which show a similar pattern of long term increasing 
activity also in the classical X-ray range: while in 2003 the 2-10 keV emission was similar 
to that seen in previous measurements with other satellites, the source flux doubled in 
the following year. These observations allowed us to study with unprecedented detail the 
properties of the persistent emission below 10 keV, by means of phase averaged and phase- 
resolved spectroscopy, and also obtain some information on the average spectral properties 
of a sample of bursts. 

2. Observations and Data Reduction 

A summary of all the XMM-Newton observations of SGR 1806-20 reported here is given 
in Table 1. We used data obtained with the EPIC instrument consisting of two MOS and 
one PN cameras (Turner et al. 2001; Struder et al. 2001). The first observation (April 
2003) was strongly affected by high particle background, resulting in a usable time of only 
~5 ks. Although SGR 1806-20 was in an apparently quiescent state, one burst could be 
identified through a careful analysis of the light curve. The next observation (October 2003) 
had a longer exposure, but no bursts were detected, despite the source was in a fairly active 
period - bursts were detected by INTEGRAL the day before and after the XMM-Newton 
observation (Mereghetti et al. 2003). 

In September 2004, SGR 1806-20 was observed again for ~50 ks while it was in an active 
state. This is the observation providing the highest quality data: no strong background flares 
were present and the source count rate was about a factor 2 higher than in the previous 
observations. Since the likely detection of bursts was anticipated, the PN camera was set 
in Small Window mode in order to reduce pile-up effects at burst peaks. In fact about 40 
bursts were detected. Finally, a Target of Opportunity observation was obtained following 
an unusual sequence of clustered soft bursts detected by INTEGRAL and Konus on October 
5 (Mereghetti et al. 2004, Golenetski et al. 2004). Again about 30 bursts were detected. 
The MOS was in Timing mode, which provides a time resolution of 1.5 ms at the expenses 
of imaging only along a single direction. 

In all the observations the medium thickness filter was used. All the data reduction was 
performed using the XMM-Newton Science Analysis Software (SAS version 6.0.0). The raw 
Observation Data Files (ODFs) were processed using standard pipeline task (epproc for PN, 
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emproc for MOS data). 

3. Spectral Analysis 

The source spectra were extracted from circular regions centered at the position of 
SGR 1806-20. For the last two observations (C and D) a radius of 40" was chosen, while for 
observations A and B, which have a higher particle background even after filtering out the 
main proton flares, a smaller radius of 25" was preferred in order to increase the signal to 
noise ratio. The background spectra were extracted from a region of the same chip as the 
source, far enough to avoid significant contamination by its photons. We selected events with 
pattern 0-4 and pattern 0-12 for the PN and the MOS, respectively. The source spectra 
were rebinned in order to have at least 30 counts per energy bin while avoid oversampling 
of the instrumental energy resolution. We limited the spectral analysis to the 1.5-12 keV 
range since, due to the high interstellar absorption (SGR 1806-20 lies in the Galactic Plane, 
at only 10° from the Galactic Center direction), very few source counts are detected at lower 
energies. 

We concentrate first on the PN spectra of the September 2004 observation (obs. C). 
This is the data set with the best statistical quality, owing to the high source count rate 
and long observing time. A fit with an absorbed power law yields a relatively high x 2 value 
(x 2 ed =1.37) and structured residuals, while a much better fit (Xr ed =0.93) can be obtained by 
adding a blackbody component. The best fit spectrum is illustrated in Fig. 1. The best fit 
parameters are photon index r=1.2, blackbody temperature kTss=0.79 keV and absorption 
Nh ~6.5x10 22 cm -2 . We have verified that the presence of bursts, which contribute to less 
than 1% of the total source counts, has a negligible effect on the source spectrum. 

In Table 2 we compare, for the four observations, the results obtained using either a 
power law or a power law plus blackbody model. Although not formally required in obs. A, 
B and D which are reasonably well fitted by single power laws, the presence of the additional 
blackbody component is consistent with all the spectra and yields systematically lower x 2 
values than the single power law fits. Except for the flux, which increased from ~1.2xl0~ n 
erg cm~ 2 s _1 of obs. A to ~2.4xl0~ n erg cur 2 s _1 of obs. C and D, all the other spectral 
parameters are, within the errors, consistent with constant values. We therefore explored 
the possibility of describing the data by forcing some of the parameters to common values 
for all the observations. This was done by jointly fitting the four data sets and resulted in 
best fit values of r=1.36, kT£#=0.65 keV and blackbody emitting area of 18 km 2 for an 
assumed distance of 15 kpc (Corbel & Eikenberry 2004). The variations can be reproduced 
by changing only the normalization of the power law component. We also analyzed in a 
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similar way the spectra obtained with the MOS cameras, obtaining results consistent with 
the PN ones. 

No evidence for spectral lines, in emission or in absorption, was found by looking at the 
residuals from the best fit models. We computed upper limits on the lines equivalent widths 
as a function of the assumed line energy and width. This was done by adding gaussian 
components to the model and computing the allowed range in their normalization. The 
most stringent results, obtained from the PN data of obs. C, are summarized in Table 3. 

4. Timing analysis and phase-resolved spectroscopy 

For the timing analysis we first corrected the time of arrival of the source events to 
the solar system barycenter and then used standard folding and phase fitting techniques to 
measure the source spin period. Pulsations were clearly detected in all the observations, 
with the period values given in Table 1. A linear fit to the four XMM-Newton period 
measurements gave a spin down F=(5.49±0.09) xlCT 10 s s _1 . 

In Fig. 2 we show the background subtracted light curves in different energy ranges for 
the four observations. Their shape, characterized by a relatively large duty cycle, is quite 
similar to that observed with RXTE (G6gu§ et al. 2002). To derive the pulsed fractions 
we fitted a sinusoid plus a constant to the background subtracted light curves. The pulsed 
fractions, defined as the amplitude of the sinusoid divided by the constant, are in the ~6-14% 
range (the values are reported in the corresponding panels of Fig. 2). 

To assess the statistical significance of possible pulse shape variations as a function 
of time and/or energy range, we compared the folded light curves using a Kolmogorov- 
Smirnov test. The results show that the difference between the soft and hard energy range 
during obs. C is highly significant (the probability that the two profiles come from the same 
underlying distribution is ~ 10~ 4 ). The corresponding value for obs. D is much higher, 0.01, 
but still indicating a possible difference between the two energy ranges. On the other hand 
the differences between the four observations are not significant (all the probabilities are 
larger than 10%). 

During obs. C enough source photons were collected to perform a meaningful phase 
resolved spectroscopy. We extracted five spectra from the phase intervals shown in Fig. 2, 
following the same method used for the average spectra. No significant spectral variations 
with phase were detected, all these spectra being consistent with the model and parameters 
of the phase-averaged spectrum, simply rescaled by an overall normalization factor. It is 
interesting to note that a fit of comparable statistical quality is obtained if the power law 
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normalization instead of the total flux is allowed to be the only phase dependent parameter. 
This suggests the possibility that the blackbody component is stable and the variations with 
the pulsation phase are due entirely to the power law component. 

Finally, the search for possible spectral lines was also extended to the phase resolved 
spectra, again with negative results. The upper limits on the equivalent widths are typically 
a factor ~3 larger than the corresponding values reported in Table 3 for the phase-averaged 
spectrum. 

5. Analysis of the bursts 

Visual inspection of the light curves binned at 0.1 s clearly showed the presence of tens 
of bursts in obs. C and D. For the other two observations the search for bursts was hampered 
by the higher background level and the worse time resolution of the PN in Full Frame mode, 
and only a single burst in obs. A could be identified. The light curves of some of the brightest 
bursts detected during obs. D are shown in Fig. 3. 

Individual bursts have too few counts for a meaningful spectral analysis. Therefore 
we extracted a cumulative spectrum of all the bursts detected during obs. C and D using 
counts from the entire PN Small Window (only the window borders have been masked). 
This corresponds to a total exposure of 12.7 s and contains about 2000 net counts in the 
2-10 keV range. The spectrum of the remaining observing time was used as background. 

Fits with simple models (power law, thermal bremsstrahlung and blackbody) all give 
formally acceptable % 2 values. The power law and the bremsstrahlung require a large ab- 
sorption (N#=10 23 cm -2 ) and the bremsstrahlung temperature is not well constrained. We 
therefore favor the blackbody model which yields an absorption value consistent with that 
of the persistent emission. The best fit parameters are kT BS =2.3±0.2 keV and N^=(6±l) 
xlO 22 cm -2 (see Fig. 4). The residuals from this best fit show a deviation at 3.3a, at 4.2 keV. 
This feature could not be reproduced in other spectra obtained with different data selections 
and binning criteria. Therefore we consider it as only a marginal evidence for an absorption 
line. Note that the absorption features previously reported in some bursts from this source 
were at a slightly different energy of ~5 keV (Ibrahim et al. 2003). 

For a few bursts included in the above analysis the large count rate at the peak can 
produce photon pile-up. An exact estimate of this effect and its correction are difficult to 
obtain in this case, owing to the rapid variability during bursts. To verify the robustness of 
our spectral results we tried different data selection aimed at reducing the fraction of piled- 
up events. We found that the shape of the spectrum does not change significantly by either 
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removing the photons of the brightest (parts of the) bursts 2 , or using an annular extraction 
region, or taking only events with pattern (single pixel). In conclusion we estimate that 
our analysis of the cumulative spectrum gives a reasonable indication of the average spectral 
shape of the bursts. 

6. Broad band spectroscopy with XMM-Newton and INTEGRAL 

The INTEGRAL observations of SGR 1806-20 (Mereghetti et al. 2005a) can in principle 
be used together with the XMM-Newton data to study the source spectrum over the broad 
1-150 keV energy range. Although the conclusions of this analysis are dependent on the 
current uncertainties in the relative calibration of the two satellites, we report here the 
results for future reference and as a possible contribution to assess this issue. 

We have considered the INTEGRAL spectra accumulated during three periods 3 which 
overlap the XMM-Newton data: 2003 March 12 - April 23, 2003 Sept. 27 - Oct. 15, and 
2004 Sept. 21 - Oct. 14. These spectra, obtained with the IBIS instrument (Ubertini et al. 
2003) and analyzed as described in Mereghetti et al. (2005a), have been fitted together with 
those of the EPIC PN. The joint fits are consistent with the parameters derived only with 
XMM-Newton (Table 2) if the IBIS data in the 20-150 keV are scaled upward by a factor 
~1.3 or ~2, respectively in the PL or PL+BB case. Alternatively, a spectral steepening 
above the EPIC range is required. For example, good fits are obtained with a photon index 
r=1.8 above 12 keV. 



7. Discussion 

Our last observation (obs. D) was carried out as a Target of Opportunity in order to 
study possible variations in the source properties related to the strong bursting activity seen 
at high energies (>20 keV) on 2004 October 5, when two clusters of powerful bursts were 
emitted in a time span of a few minutes (Mereghetti et al. 2004). The individual events had 
spectra and duration similar to those of the normal bursts, but they had rather high peak 
fluxes and were grouped in a way never observed before in this source, yielding a total fluence 



2 the best fit blackbody model excluding data with more than 20 counts/frame gives kTss=2.32±0.23 
keV and N H =(6.0±1.3) x 10 22 cni~ 2 

3 since long exposure times are required to detect the source in the hard X-ray range, it is not possible 
to restrict the analysis to the INTEGRAL data strictly simultaneous with the XMM-Newton observations 
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of ~10~ 4 erg cm -2 (Golenetskii et al. 2004). The XMM-Newton observation, which started 
only 27 hours after this energetic event, did not show particularly striking changes in the 
source spectrum and pulse profile. There is some evidence of a more structured pulse profile 
above 5 keV, but the statistical significance is not compelling. This means that the October 
5 event did not cause changes in the mechanisms responsible for the persistent emission and 
in magnetic field configuration that were strong and/or long lasting enough to affect the 
source properties significantly after ~1 day. 

Thanks to the high statistics in the spectrum of the persistent emission obtained in 
September 2004 (obs. C), we could set stringent limits on the presence of lines and found 
evidence for an additional blackbody component with temperature kT## ~0.8 keV. This 
component was not required in the spectra from the previous observations, when the source 
had a lower luminosity, but the data are consistent with it being always present with a 
constant temperature and emitting area. A similar thermal component, although with a 
lower temperature was also reported in Chandra observations of SGR 1806-20 performed in 
May- June 2004 (Woods et al. 2004). During these observations the source soft X-ray flux 
was already larger than the "historical" level of ~1.3xl0 -11 erg cm~ 2 s _1 . A blackbody 
component has also been observed in the other well studied soft repeater SGR 1900+14 
(Woods et al. 2001), where it was more prominent when the source was in a relatively 
quiescent state. The two components spectrum reinforces the similarity between SGRs and 
Anomalous X-ray Pulsars, where the two component model has an ubiquitous character (e.g. 
Mereghetti et al. 2002a), and which are also believed to be magnetars. 

The power law photon indeces 4 of SGR 1806-20 obtained in the four XMM-Newton 
observations (r ~1.5) are significantly smaller than those observed in 1993 with ASCA 
(r=2.2±0.2; Sonobe et al. 1994) and in 1998-2001 with BeppoSAX (r=1.97±0.09; Mereghetti 
et al. 2002b). This indicates that a spectral hardening occurred between September 2001 
and April 2003. In this time interval, only RXTE observed SGR 1806-20 during a monitoring 
program mainly focussed on the source timing properties and no spectral results (which could 
establish when exactly the spectral change occurred) have been reported to date. On the 
other hand, these RXTE observations indicate that the average spin-down rate changed in 
2000. While the early sparse period measurements with ASCA and BeppoSAX (Mereghetti 
et al. 2002b), as well as a phase-connected RXTE timing solution spanning February-August 
1999 (Woods et al. 2000), are consistent with an average P ~ 8.5 x 10~ n s s _1 , subsequent 
RXTE data indicate a spin-down larger by a factor ~4 (Woods et al. 2002). Our period 
measurements show a further increase in the average P, as shown in Fig. 5. 



4 since the previous results were described without the additional blackbody component, we compare the 
fits with a single power law 
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A correlation between spectral hardness and spin-down rate is present in the sample of 
AXPs and SGRs. The sources with the harder spectrum have a larger long term spin-down 
rate (Marsden & White 2001). The results shown in Fig. 5 for SGR 1806-20 indicate, for 
the first time, that such a correlation also holds within different states of a single source. 

All the results discussed above are consistent with the magnet ar scenario. In particular 
the observed long term variations, the correlation between spectral hardening and spin-down 
rate, and the increase in the bursting activity fit reasonably well with what is expected in a 
twisted magnetosphere (see Thompson, Lyutikov & Kulkarni 2002). In this scenario, SGRs 
and AXPs differ from standard radio pulsars since their magnetic field is globally twisted 
inside the star, up to a strength about 10 times the external dipole, and, at intervals, can 
twist up the external field. Twisted, force-free magnetospheres (B^ ^ 0) support currents 
and the charges flowing in the magnetosphere may produce a significant optical depth to 
electron and ion scattering. Since scattering is resonant at the cyclotron frequency u>b, 
thermal photons emitted at the star surface scatter at different radii, where ui = ujb{R)- In 
the case of electrons, the charge distribution is spatially extended and repeated scatterings 
lead to the formation of a high-energy tail (instead of a narrow line). A gradually increasing 
twist results in a larger optical depth and this causes a hardening of the X-ray spectrum. 

At the same time, the spin down rate increases because, for a fixed dipole field, the 
fraction of field lines that open out across the speed of light cylinder grows. Since both 
the spectral hardening and the spin-down rate increase with the twist, the model predicts 
that they should be correlated. This trend is actually present in the data reported here (see 
Fig. 5). According to the magnetar model, the stresses building up in the neutron star crust 
and the magnetic footprints movements can lead to crustal fractures which can be energetic 
enough to explain the observed increase in the bursting activity (but see Jones 2003). 

The charges present in the magnetosphere also provide a large optical depth to resonant 
proton (ions) cyclotron scattering. The proton resonance sits much closer to the star surface 
than the electron resonance, and therefore is less sensitive to the broadening caused by the 
radial dependence of B. However, this is not sufficient to make spectral lines detectable 
because, in the general case, positive charges are not confined in a thin layer close to the star 
surface. This implies that lines in the persistent emission are difficult to observe because the 
scattering occurs in a region where the magnetic field is varying. 

Interestingly, this scenario may provide an explanation for the transient appearance of a 
cyclotron line as that tentatively detected at 4 keV during the bursts seen with XMM-Newton 
and at 5 keV with RXTE (Ibrahim et al. 2002; 2003), at times in which the source has a 
particularly large luminosity and hard spectrum. In fact, when the luminosity at ub(R*) 
exceeds the luminosity produced by surface heating by returning currents (L r x ~ 10 35 - 
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10 36 erg s -1 , see eq. [34] of Thompson, Lyutikov & Kulkarni 2002) protons may accumulate 
close to the surface owing to the large radiation force. Under these conditions, i.e. the 
existence of a twisted, current-carrying magnetosphere and of a large transient flux (above 
L r x) at the cyclotron resonance, the formation of a line may be expected. The persistent 
luminosity of SGR 1806-20 in the range 4-5 keV is below L r x, while the average burst 
luminosity in the same energy range is > 2 x 10 36 erg s _1 , quite larger than U^. The different 
energy at which the line has been possibly detected with XMM-Newton and RXTE may be 
explained in terms of a variation 10%) of the average surface field along the line of sight 
in the two epoches, which in turn can be due to different twisting properties or quadrupolar 
components, to which the proton cyclotron line is very sensitive. 

Up to now limited spectral information has been obtained for SGR bursts below 20 
keV, in particular with high spectral resolution and good sensitivity at a few keV. Several 
studies have provided increasing evidence that the optically thin thermal bremsstrahlung 
model, which gives a good phenomenological description of the SGR burst spectra in the 
hard X-ray range, is inconsistent with the data below 15 keV (Fenimore, Laros & Ulmer 
1994; Olive et al. 2004). Recently Feroci et al. (2004) analyzed a sample of bursts from 
SGR 1900+14 using BeppoSAX data in the 1.5-200 keV range and obtained good fits with 
a model consisting of two blackbodies with average temperatures of about 3 and 10 keV. 
Our results are consistent with a similar situation also for SGR 1806-20 and confirm earlier 
indications in this sense (e.g., Laros et al. 1986). 



8. Conclusions 

The XMM- Newton monitoring of SGR 1806-20 carried out in the last two years has 
allowed us to perform a detailed study of its soft X-ray emission, in particular for what 
concerns the spectral properties and long-term variability. Our results are particularly in- 
teresting when compared to earlier ASCA and BeppoSAX measurements, and put in the 
broader context of observations at higher energies. The emerging scenario is that of a long- 
term growth in the level of non-thermal magnetospherical activity, finally leading to the 
exceptional event observed on 2004 December 27. 

More specifically, we found that: 

• a single power law is inadequate to fit the 2-10 keV spectrum with the highest statistical 
quality: an additional blackbody component, too weak to be detected in less sensitive 
observations, is required 

• the temperature of such a component (kT BB ~ 0.8 keV) is slightly higher than that 
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observed with Chandra three months earlier (kT B B ~ 0.54 keV, Woods et al. 2005) 

• although not formally required in the fits, a blackbody component with constant tem- 
perature and luminosity is compatible with all the XMM-Newton observations; this 
requires that the variations in the source luminosity and spectrum, both on long term 
and as a function of the pulsar phase, be accounted for by changes in the power law 
component 

• no absorption or emission lines are present in the 2-10 keV persistent emission with 3a 
upper limits from 30 to 100 eV on the equivalent width 

• the 2-10 keV luminosity nearly doubled during the first half of 2004, reaching a level 
of ~ 10 36 erg s -1 ; the same trend has been observed above 20 keV with INTEGRAL 

• the spectrum has been monotonically hardening during the last decade, with the power 
law photon index decreasing from 2.2 to 1.5 

• the increase in the spectral hardness and in the average spin-down rate follow, within 
the same source, the correlation previously found by comparing different SGRs and 
Anomalous X-ray Pulsars 

• these changes were not accompanied by large variations in the pulse profile, even after 
the relatively strong sequence of bursts of 2004 October 5 

• a blackbody with temperature of about 2 keV provides a good fit to the average burst 
spectra below 10 keV 

• there is marginal evidence for an absorption feature at 4.2 keV in the bursts spectra 
which needs to be confirmed with further data 



These findings fit reasonably well in the magnetar scenario in which the effects of a 
twisted magnetosphere are considered (Thompson et al. 2002). In this context we have also 
proposed an interpretation for the possible lines observed during the bursts. 

It is clear that the hard X-ray bursting activity and the properties of the so called 
"persistent" emission, until recently only observed below 20 keV, are strongly connected. 
Dramatic changes in the source properties are expected as a consequence of the enormous 
energy release occurred during the 2004 December 27 event. The total isotropic energy in 
the hard X-ray range emitted in this giant flare was of a few 10 46 erg (Mazets et al. 2005; 
Terasawa et al. 2005; Hurley et al. 2005; Mereghetti et al. 2005b). A comparison of future 
X-ray observations of SGR 1806-20 with the results reported here for the "pre-Giant-Flare 
epoch" will certainly lead to a better understanding of this class of sources. 
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Table 1: Journal of XMM-Newton observations and timing results. 

Obs. Date Duration Mode (a) and exp. time and exp. time Pulse Period 

(ks) PN camera MOS cameras (s) 

A 2003 Apr 03 30 FF (5.4 ks) LW (6 ks) 7.5311±0.0003 

B 2003 Oct 07 22 FF (13.4 ks) LW (17 ks) 7.5400±0.0003 

C 2004 Sep 06 52 SW (36.0 ks) LW (51 ks) 7.55592±0.00005 

D 2004 Oct 06 19 SW (12.9 ks) Ti (18 ks) 7.5570±0.0003 



a FF = Full Frame (time resolution 73 ms); LW = Large Window (time resolution 0.9 s); SW = Small Window 
(time resolution 6 ms); Ti = Timing (time resolution 1.5 ms) 



Table 2. Summary of the spectral results^ 
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1.01 
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(a) Errors are a ^ t ne go% c.l. for a single interesting parameter 
( 6 ) Radius at infinity assuming a distance of 15 kpc 
( c ) Absorbed flux of the power law component as observed by the PN in the 2-10 keV energy range 
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Table 3: Upper limits (at 3 a) on spectral features in the PN spectra of SGR 1806-20 from 
observation C. 



Energy range 


a (eV) 


Equivalent width 


2-3 keV 





< 28 eV 




100 


< 65 eV 




200 


< 100 eV 


3-7 keV 





< 22 eV 




100 


< 32 eV 




200 


< 55 eV 


7-10 keV 





< 55 eV 




100 


< 70 eV 




200 


< 110 eV 
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Fig. I.— EPIC PN spectrum of SGR 1806-20 for the observation of 2004 September 6 (obs. 
C). Top: data and best fit model (solid line). The two spectral components, power law and 
blackbody (lower curve) are indicated by the dashed lines. Bottom: residuals from the best 
fit model in units of a. 
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Fig. 2.— Folded light curves in the total (2-10 keV), soft (2-5 keV), and hard (5-10 keV) 
energy range for the four observations. The background has been subtracted. The pulsed 
fraction is indicated on each panel (la errors). The vertical lines in obs. C indicate the phase 
intervals used for the phase resolved spectroscopy. 
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Fig. 3. — Light curves of four bursts detected during obs. D. Less than 0.02 background 
events per frame are expected. 
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Fig. 4. — EPIC PN spectrum of the sum of 69 bursts detected in obs. C and D. Top: data 
and best fit model with a blackbody of temperature kT=2.3 keV. Bottom: residuals from 
the best fit model in units of a. 
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Fig. 5. — Long term evolution of the pulse period (top) and power law photon index (bottom) 
of SGR 1806-20 . To derive the fourth BeppoSAX data point in the lower panel we have 
analyzed an unpublished observation performed on 6 September 2001. 



